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摘 要 : 黄土 高 原 是 中 国 重 要 的 生态 屏障 。 在 全 球 变 化 的 背景 下 ,黄土 高 原 的 碳 收 支 平衡 备 受 关注 。 本 研究 基于 


MODI7A3HGF 数据, 通过 GSMSR 模型 趋势 分 析 、 差 异 分 析 以 及 地 理 探测 器 等 方法 ,分 析 2000 一 2020 年 黄土 高 原 
的 碳 源 / 汇 特征 ,揭示 2000 一 2020 年 该 区 域 生 态 系 统 净 生产 力 (NEP) 的 时 空格 局 及 其 驱动 因素 。 同 时 ,将 研究 区 按 
经 度 方向 划分 为 西 . 中 、 东 三 个 子 区 域 ,比较 不 同 区 域内 驱动 因素 的 差异 性 。 结 果 表 明 :(1) 近 20 a 黄土 高 原 
49.69% 的 区 域 从 碳 源 向 碳 汇 转变 ; NEP 随时 间 波 动 上 升 ,在 东南 部 高 于 西北 部 ,多 年 平均 值 为 12g C.m?.a ^ (2) 水 


分 条 件 是 影响 NEP 空间 分 布 的 主要 自然 因素 ,而 土地 利 ) 


] 类 型 则 是 影响 NEP 空间 分 布 的 主要 人 为 因素 ;不 同 因子 


间 的 交互 作用 对 NEP 的 影响 普遍 大 于 单个 因子 ,(3) 西 .中 、 东 三 个 子 区 域 NEP 的 驱动 因子 存在 明显 的 空间 分 异 特 


响 ,其 中 ,以 土地 利用 类 型 为 代表 的 人 为 干扰 最 强 。 


征 ,中 、 西 部 地 区 受气 候 影 响 较 多 ,以 降水 .湿度 等 水 分 条 件 为 主 ;东部 地 区 受 地 形 .气候 .人 类 活动 等 因素 综合 影 
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植被 是 陆地 生态 系统 的 重要 组 分 之 一 ,在 调节 
碳 平衡 方面 起 着 至 关 重 要 的 作用 ,增强 植被 左 汇 
功能 是 减缓 大 气 CO; 浓度 上 升 和 全 球 变 暧 的 有 效 途 
径 , 是 实现 “ 碳 中 和 "目标 的 关键 因素 ”。 生 态 系统 
净 生 产 力 (NEP) 是 指 陆地 生态 系统 与 大 气 的 净 碳 交 
换 ” ,其 描述 了 一 个 生态 系统 在 单位 时 间 内 能 够 固 
定 大 气 中 CO 的 量 , 是 评估 生态 系统 碳 源 / 汇 功能 的 
重要 指标 “。 已 有 人 研究 表明 ,陆地 生态 系统 具有 明 
显 的 矶 汇 特征 ,但 是 , 碳 汇 的 空间 分 布 存 在 较 大 差 
异 , 且 易 随 外 部 条 件 的 改变 而 改变 。 区 域 生 态 系 统 
的 碳 源 / 汇 变化 及 其 驱动 机 制 的 确定 仍 具 有 挑战 性 ， 
因此 ,区 域 NEP 佑 算 、 时 空格 局 及 其 驱动 因子 分 析 
也 成 为 当前 全 球 变化 研究 的 热点 问题 一 。 

黄土 高 原 具 有 地 形 复杂 气候 多 变 等 突出 特 
点 ,在 我 国生 态 环境 保护 和 社会 经 济 发 展 中 都 发 挥 
着 重要 的 作用 。 长 期 以 来 该 地 区 由 于 人 类 活动 频 
FE ,植被 退化 严重 ,水 土 流失 剧烈 发 生 ”" ,是 我 国 
十 分 典型 的 生态 环境 脆弱 区 。 近 20 a, 由 于 “退耕 还 
林 ( 草 )" 工 程 " 的 实施 ,该 地 区 植被 覆盖 有 所 增加 ， 
生态 环境 显著 改善 。 目 前 ,有 关 黄 土 高 原 归 一 化 植 
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PAGEL (NDVI) 植被 覆盖 度 (FVC) 以 及 植被 净 
初级 生产 力 (NPP)™ 等 指标 进行 区 域 生态 环境 监测 
和 评估 的 研究 较 多 ,但 进一步 探究 植被 碳 源 / 汇 时 空 
分 布 的 研究 相对 较 少 ,而 有 关 碳 源 / 汇 时 空 分 布 驱动 
因子 的 研究 更 为 少见 , 既 有 的 研究 中 也 存在 如 下 问 
题 ; 一 是 驱动 因子 的 选取 相对 较 少 ,对 NEP 的 空 
间 格 局 的 解释 力 有 限 ;二 是 研究 仅仅 分 析 了 驱动 因 
子 与 NEP 的 相关 性 ,未 能 量化 驱动 因子 对 NEP 的 贡 
献 ; 三 是 局 限于 单 因子 对 NEP 的 贡献 ,忽视 了 不 同 
子 之 间 对 NEP 贡献 的 协同 作用 ;四 是 过 于 强调 黄 
土 高 原 的 整体 性 ,而 针对 内 部 的 差异 性 的 研究 相对 
不 足 ; 咪 需 开 展 更 为 深入 的 人 研究。 

针对 上 述 问 题 ,本 研究 以 大 的 省 级 行政 区 为 划 
分 单元 ,将 研究 区 自 西向 东 划 分 为 三 个 子 区 域 ,以 
探究 研究 区 内 部 的 差异 性 。 通 过 NEP 估算 模 型 对 
研究 区 2000—2020 年 NEP 进行 估算 。 采 用 Theil- 
Sen Median #4344} Wr , Mann- Kendall 显著 性 检验 以 
及 差异 分 析 探 究 研 究 区 NEP 的 时 空 变 化 。 最 后 利 
用 地 理 探测 器 对 人 研究 区 内 影响 NEP 的 驱动 因子 进 
行 筛选 和 分 析 ,并 定量 评估 其 对 NEP 的 影响 水 平 以 
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及 不 同 驱动 因子 之 间 的 交互 作用 。 本 研究 由 在 探 ”平台 中 重 采样 为 1 km。 


究 黄 土 高 原 NEP 的 时 空 演变 特征 ,揭示 2000 一 2020 
年 研究 区 生态 建设 改善 的 效果 ,为 该 区 域 可 持续 发 
展 、 生 态 环境 保护 提供 科学 依据 和 决策 参考 。 


1 研究 区 概况 与 方法 


1.1 研究 区 概况 

黄土 高 原 位 于 黄河 中 上 游 (33"43'~41?16'N， 
100*54' -114?33' E) ,地 处 由 沿海 向 内 陆 , 平 原 向 高 
原 过 渡 地 带 。 年 平均 降水 量 从 西北 向 东南 增加 
(200-750 mm) ,集中 在 7 一 9 月 ,多 以 暴雨 形式 出 
现 。 莹 发 量 一 般 高 于 降水 量 。 年 平均 气温 为 4.3~ 
14.3 C, DXA HOSEA PY te RA) AX Er o 1) AR LI 
地 、 黄 土 沟 密 和 河谷 平原 。 研 究 区 按 经 度 方向 自 西 
向 东 划 分 为 西 .中 、 东 三 部 分 ,西部 地 区 地 处 干旱 半 
干旱 区 ,涵盖 了 宁夏 全 域 及 青海 和 甘肃 的 部 分 地 
区 ,区 内 海拔 多 在 2000~3000 m, 是 黄土 高 原 地 热 最 
高 的 地 区 ;中 部 地 区 主体 为 关中 平原 和 陕 北 高 原 ， 
降水 分 布 具 有 纬度 地 带 性 ,区 南 为 半 湿 润 区 ,区 北 
为 半 干 旱 区 ;东部 地 区 覆盖 了 山西 全 域 和 河南 的 部 
分 地 区 ,地 势 东北 高 西南 低 ,地 形 地 貌 复 杂 , 多 为 山 
地 丘陵 ,气候 干 湿 状 况 主 要 为 半 湿 润 区 ,是 黄土 高 
原 降 水 最 充沛 的 区 域 (图 1)。 
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图 1 研究 区 范围 及 高 程 
Fig. 1 Range and elevation ofthe study area 


1.2 数据 来 源 

本 研究 收集 了 研究 区 的 NPP 和 环境 因子 数据 
( 表 1)。 目 前 ,NPP 数 据 集 已 广泛 应 用 于 估算 碳 源 / 
汇 ”。 空 间 数 据 的 空间 参考 设置 为 CGCS2000_ 
3_Degree_GK_CM_105E ,空间 分 辩 率 在 ArcGIS 10.8 


1.3 研究 方法 

1.3.1 NEP 估 算 模 型 NEP 指 的 是 生态 系统 碳 循环 
的 净 吸 收 量 或 净 储 存量 ,已 被 广泛 用 于 估算 陆地 生 
态 系统 碳 源 / 汇 。 若 NEP>0 ,表示 植被 固定 的 碳 量 高 
于 土壤 微生物 呼吸 排放 的 碳 量 , 生 态 系统 表现 为 碳 
汇 ; 若 NEP<0 ,表示 植被 固定 的 碳 量 低 于 土壤 微 生 
物 呼 吸 排放 的 碳 量 , 生 态 系统 表现 为 碳 源 "”"。 年 
NEP 被 定义 为 年 NPP 和 年 土壤 异 养 呼吸 (Rh) 的 差 
,其 中 Rh 是 由 Rh 和 土壤 自 养 呼 吸 (Rs) 之 间 关 系 
的 拟 合 方程 来 确定 。Rh 对 降水 温度 和 土壤 有 机 碳 
(SOC) 敏 感 ,并 且 与 局 成 比例 。 根 据 1994 一 2013 年 


的 多 年 记录 ,Chuai 等 "计算 出 适合 中 国 地 区 的 Rh 
和 Rs 的 线性 拟 合 方程 。 

NEP = NPP - Rh (1) 

Rh =0.4679 x Rs + 114.42 (2) 


根据 Yu 等 2 构建 的 GSMSR 模型 ,研究 区 月 土 


式 确定 : 
RS month = (0.588 + 0.1 18 x SOC)e 
x(P 2.972) + (P+ 5.657) x 30 


n(1.83x e" *7)x T+ 10 


(3) 
12 

Rs annual = ` Rs month ( 4 ) 
i=l 


IRH : SOC 为 0~20 cm 土壤 有 机 碳 密度 (kg Cm?) ;P 
为 月 降水 量 (cm );7 为 月 平均 温度 (% )。GSMSR 已 
被 证 明 在 模拟 中 国 范围 内 的 Rs 是 可 靠 的 ,因此 本 文 
采用 该 模型 模拟 2000 一 2020 年 黄土 高 原 的 NEP。 
1.3.2 趋势 分 析 Theil-Sen 趋势 中 位 数 分 析 是 一 种 
稳健 的 非 参数 检验 的 趋势 分 析 方法 2 ,计算 效率 
高 ,并 且 已 经 在 之 前 的 研究 中 得 到 广泛 应 用 ”。 该 
表达 式 可 写 为 : 


NEP, - NEP, 


B= Medion pnr 
式 中 : B 表示 NEP 的 变化 趋势 ;NEP 和 NEP 分 别 表 
示 年 份 : 和 7 的 时 间 序 列 NEP 值 (本 研究 中 2000< i< 
j $2020) ,表示 研究 周期 的 长 度 。 如 果 B>0 表示 
NFEP 呈 上 升 趋势 ,而 8<0 则 表示 NEP 呈 下 降 趋 势 。 
Mann-Kendall 趋势 检验 是 一 种 基于 等 级 的 非 参 
数 检验 ,可 以 用 来 检验 NEP 变 化 趋势 的 显著 性 的。 
具有 不 需要 样本 服从 归 一 化 分 布 并 且 不 受 异 常 值 
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表 1 数据 来 源 
Tab. 1 Data sources 
数据 时 间 A BE 数据 来 源 
净 初 级 生产 力 (NPP) 2000 一 2020 年 500m 美国 国家 航空 航天 局 (https:Wladsweb.modaps.eosdis.nasa.gov ) 
海拔 (ALT) B 30m 地 理 空间 数据 云 (http://www.gscloud.cn ) 
年 蒸发 量 (EVP) 2000 年 ,2010 年 ,2020 年 1 km 
年 平均 地 温 (GST) 2000 年 ,2010 年 ,2020 年 1 km 
年 平均 气压 (PRS) 2000 年 ,2010 年 ,2020 年 1 km 
年 平均 相对 湿度 (RHU) 2000 年 ,2010 年 ,2020 年 1 km 
年 日 照 时 数 (SSD) 2000 年 ,2010 年 ,2020 年 1 km 中 国 科学 院 资源 环境 科学 与 数据 中 心 (https://www.resdc.cn) 
国内 生产 总 值 (GDP) 2000 年 ,2010 年 ,2019 年 1 km 
人 口 密度 (PD) 2000 年 ,2010 年 ,2019 年 1 km 
年 平均 风速 (WIN) 2000 年 ,2010 年 ,2020 年 1 km 
土壤 类 型 (SOT) 1995 年 1 km 
E 量 年 
ire ) a 在 . 国家 地 球 系统 科学 数据 中 心 (https://www.geodata.cn) 
土地 利用 / 窗 盖 (LULC) 2020 年 30m Zenodo(https://zenodo.org ) 
坡度 (SLO 一 30 m DEM 坡 度 分 析 
坡 向 (ALP) 一 30 m DEM 坡 向 分 析 
与 不 透水 面 距离 (DFT) 2020 年 1 km LULC 提 取 
与 道路 距离 (DFR ) 2020 年 1 km LULC 提 取 
与 水 系 距 离 (DFW) 2020 年 1 km LULC 提 取 
0-20 cm 土壤 碳 密度 (SOC) ”1995 年 = 国家 青藏 高 原 科 学 数据 中 心 (http:Wdata.tpdc.ac.cn ) 


干扰 的 优点 。 
1.3.3 差异 分 析 通过 逐 像素 计 算 2000 年 和 2020 
年 NEP 的 差 值 ,得 到 了 黄土 高 原 NEP 的 变化 幅度 ，， 
D, - NEP; - NEP; (6) 

式 中 :D; 为 2020 年 和 2000 年 第 ; 行 第 j 列 像素 的 NEP 
之 差 ;4 为 2020 年 ;为 2000 年 ; NEP; 为 2020 年 第 i 
行 第 j 列 像素 点 的 NEP 值 ; NEP; 为 第 i 行 第 j 列 像素 
点 在 2000 年 的 NEP 值 。 
13.4 地 理 探测 器 ”地理 探测 器 是 一 种 空间 统计 方 
法 ,基于 四 个 模块 揭示 空间 异 质 性 并 量化 驱动 因素 
对 响应 变量 的 影响 “”。 本 文采 用 因子 探测 顺和 交 
互 探测 器 分 析 研 究 区 以 及 其 内 部 区 域 NEP 的 驱动 
因子 。 

(1) 因子 探测 器 ,NEP 的 驱动 因素 的 解释 力 可 
以 通过 因子 探测 需 的 4 值 来 衡量 。4 值 越 大 ,对 NEP 
空间 分 布 的 解释 力 越 强 。 表 达 式 如 下 : 


(7) 


式 中 :hn=1,2,3,… 沁 表示 相关 或 独立 因素 的 层次 ; 
NN 和 NN 分 别 表示 层 内 单元 数 和 整个 区 域 中 的 单元 


数 ; o; 和 o* 分 别 是 hh 层 的 方差 和 整个 区 域 的 Y 值 的 
方差 。 

(2) 交互 探测 器 ,交互 探测 器 可 以 通过 比较 两 
^i vr PER XS [a OG ] X3 [g Xa) EJ X IL XS HY 8. 
作用 [g(XNK)] 的 g 统 计 值 ,来 确定 两 个 驱动 因素 的 
解释 力 是 增强 ,减弱 还 是 相互 独立 。 交 互 类 型 的 详 
细 信 息 见 表 2。 


表 2 驱动 因子 的 交互 类 型 


Tab.2 Interaction types of driving factors 


交互 类 型 交互 关系 
非 线 性 减弱 qGOXX;)«Min[q CX.) ,4(X;)] 
单 因子 非 线 性 减弱 ”Min[g(Z) gX) qX X) «Max[q X) , 
q(Q6)] 
双 因 子 增强 q GG)» Max[q X) ,gq(X)] 
独立 qGDX:)2q 06 +g) 
非 线性 增强 q(XiNX2)>q( Xi) *q (X2) 


2 结果 与 分 析 


2.1 NEP 时 间 变 化 特征 

2000—2020 年 黄土 高 原 NEP 多 年 平均 值 为 
12.00 gC.m?*…a'。 最 小 值 出 现在 2001 年 ,为 -87.11 
g Cem? a” OCA HUE 2018 4E , 73 88.80 g Cem? 
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a (图 2d)。 在 内 部 三 个 区 域 中 ,中 部 NEP 多 年 平均 
值 最 高 ,为 44.91 g Cema”, PIXE NEP LEFFE 
最 小 , 仅 为 -7.46 g Cm a, KAB NEP 多 年 平均 值 
介 于 二 者 之 间 , 为 10.09 g C-m?-a, 研究 区 NEP 年 
际 变化 特征 明显 ,变化 速率 为 7.07 gC.m*…a!。 其 
中 ,中 部 地 区 NEP 上 升 速率 最 高 ,达到 了 8.99 g C- 
nm2.a。 西 部 和 东部 地 区 NEP 也 呈 上 升 趋势 ,变化 
速率 分 别 为 6.17 gC.m2.a 和 6.69 g Cem”. a" Cd 
2a~ 图 2c)。 
2.2 基于 NEP 的 碳 源 / 汇 特征 分 析 

2000 年 黄土 高 原生 态 系 统 具 备 碳 源 / 汇 特征 的 
区 域 面积 占 比分 别 为 78.80% 和 21.20%。 具 备 碳 源 
特征 的 区 域 主要 分 布 在 西北 地 区 以 及 关中 平原 、 汾 
河谷 地 。 具 备 碳 汇 特征 的 区 域 主要 分 布 在 吕梁 山 、 
秦岭 (图 3a)。2020 年 黄土 高 原 具备 碳 源 / 汇 特征 的 
区 域 面积 占 比分 别 为 29.20% 和 70.80%。 具 备 碳 源 
特征 的 区 域 面积 大 幅 减 少 , 主 要 分 布 在 研究 区 西北 
缘 及 各 大 中 城市 (图 3b)。 总 体 上 呈现 出 碳 源 西北 ， 
碳 汇 东南 的 格局 。 从 2000—2020 年 黄土 高 原 有 


75 p (a) 西部 地 区 

50L 
个 25+ 
a 
g 0 | 
O 
2 25} 
A 
Z sol 

aL y=6.17x-75.27 
R?-0.85 
-100 
2000 2005 2010 2015 2020 


90 r- (c) 东部 地 区 


49.78% 的 区 域 发 生 了 碳 源 和 碳 汇 的 转换 。 其 中 ,从 
具备 碳 源 特 征 的 区 域 转换 为 具备 碳 汇 特征 的 区 域 
面积 占 比 为 49.69% ,从 碳 汇 特征 的 区 域 转换 为 碳 源 
特征 的 区 域 面 积 占 比 仅 为 0.09%( 图 3c, 图 3d)。 
2.3 NEP 趋 势 分 析 与 差异 分 析 

黄土 高 原 85.01% 的 区 域 NEP 呈 上 升 趋势 , 极 显 
著 上 升 区域 占 55.17% ,广泛 分 布 于 研究 区 。 黄 土 高 
原 0.96% 的 NEP 呈 下 降 趋 势 ,分 布 在 各 大 中 城市 (图 
4b)。 在 人 研究 区 内 部 ,中 部 地 区 NEP 星 上 升 趋 热 的 
区 域 占 比 最 大 ,为 94.83%, 极 显著 上 升 区 域 占 中 部 
地 区 的 80.67% ,是 研究 区 NEP 上升 最 为 明显 的 区 域 
( 表 3)。 西 部 地 区 NEP 旺 上 升 趋势 的 区 域 占 西部 地 
区 的 86.07% , 极 显著 上 升 的 区 域 占 比 为 50.18% , 仅 
次 于 中 部 地 区 。 东 部 地 区 NEP 呈 上 升 趋势 的 区 域 
占 东 部 地 区 的 76.65% , 呈 极 显著 上 升 的 区 域 占 比 为 
42.24%。 相 较 于 中 .西部 地 区 而 言 ,东部 地 区 NEP 
上 升 趋势 相对 较 缓 ,以 相对 稳定 的 区 域 居多 , 占 比 
为 22.46%。 

黄土 高 原 2020 年 NEP 大 于 2000 年 NEP 的 像 元 
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图 2 2000 一 2020 年 黄土 高 原 及 其 子 区 域 NEP 的 年 际 变化 
Fig. 2 Annual variation of NEP in Loess Plateau and its sub-regions from 2000 to 2020 
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(b) 2020 年 NEP 空 间 格局 
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(d) 2000 一 2020 年 黄土 高 原 碳 源 和 碳 汇 面积 变化 
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图 3 2000 一 2020 年 黄土 高 原 NEP 空 间 格局 及 碳 源 汇 转化 
Fig.3 NEP spatial pattern and carbon source-sink transformation in Loess Plateau from 2000 to 2020 
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(b) 2000 一 2020 年 NEP 趋 势 分 析 
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图 4 2000 一 2020 年 NEP 差 异 分 析 和 趋势 分 析 
Fig. 4 Difference analysis and trend analysis of NEP from 2000 to 2020 


占 比 为 96.79%( 表 4) , NEP 差 值 小 于 0 的 区 域 主要 
分 布 在 各 省 会 或 其 他 中 心 城市 (图 4a)。NEP 差 值 
介 于 0~300 的 区 域 占 比 以 西部 地 区 居多 ,为 
98.44%。 中 部 地 区 NEP 差 值 大 于 300 的 区 域 占 比 相 


对 最 大 ,为 4.33%。 此 外 ,NEP 差 异 最 明显 的 区 域 呈 
现 出 从 甘肃 和 陕西 的 交界 处 向 周围 扩散 的 分 布 
特征 。 

综 上 ,人 研究 时 段 内 ,黄土 高 原 NEP 均 以 上 升 态 
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表 3 2000 一 2020 年 NEP 变 化 趋势 的 统计 
Tab.3 Statistics for 2000-2020 trend analysis of NEP 


斜率 IZI 值 NEP 显 著 性 水 平 FERME 
西部 地 区 中 部 地 区 东部 地 区 研究 区 
<0 22.58 极 显著 减少 0.01 0.01 0.03 0.02 
«0 =1.96 显著 减少 1.44 0.26 0.86 0.94 
- «1.96 相对 稳定 12.47 4.91 22.46 14.03 
>0 三 1.96 显著 增加 35.89 14.16 34.41 29.84 
>0 =2.58 极 显 著 增加 50.18 80.67 42.24 55.17 


表 4 2000 一 2020 年 NEP 差异 分 析 的 统计 
Tab. 4 Statistics for 2000-2020 difference analysis of NEP 


NEP 差异 水 平 百分比 /% 
KgC'm^a?) 西部 地 区 中 部 地 区 ”东部 地 区 ”研究 区 
<0 0.81 1.82 3.70 3.21 
0-100 86.51 1347 22.65 26.74 
100-200 6.22 40.83 61.13 42.88 
200-300 5.71 39.85 12.14 24.60 
>300 0.75 4.33 0.38 2.57 


势 为 主 , 中 部 地 区 NEP 呈现 上 升 态 势 的 区 域 占 比 和 
增加 程度 均 优 于 其 他 地 区 ,西部 地 区 次 之 ;东部 地 
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[X NEP 呈 相 对 稳定 的 区 域 占 比 (22.46% ) 大 于 中 
(4.9196) PE (12.47% ) 部 地 区 ;NEP 下 降 的 区 域 主要 
分 布 在 大 中 城市 或 其 他 经 济 发 达 区 域 。 
2.4 NEP 的 驱动 因子 
2.4.1 因子 探测 ”黄土 高 原 及 子 区 域 NEP 和 18 个 
潜在 驱动 因子 经 因子 探测 器 得 到 了 初步 结果 。 为 
了 便于 分 析 表 达 , 本 研究 将 不 同 区 域 的 潜在 驱动 因 
子 按 g 值 进行 排序 ,选取 累计 贡献 值 超过 80% 的 潜 
在 驱动 因子 进行 后 续 的 交互 作用 探测 。 

通过 因子 探测 器 研究 发 现 (图 $) ,降水 因子 对 
人 研究 区 内 NEP 空 间 分 布 的 解释 力 最 强 ,贡献 最 大 


(b) 中 部 地 区 


(d) 黄土 高 原 


图 5 黄土 高 原 及 其 子 区 域 NEP 因子 探测 结 


Fig. 5 NEP factor detection results of Loess Plateau and its sub-regions 
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(q-0.4) ,相对 湿度 的 贡献 次 之 (q4=0.31)。 蔡 发 量 、 
土地 利用 类 型 .土壤 类 型 及 日 照 时 数 对 NEP 的 贡献 
位 于 中 等 水 平 (g >0.2)。 在 空间 上 ,驱动 因子 对 区 
域 NEP 的 贡献 呈现 出 明显 的 空间 分 异 格局 。 降 水 ， 
蒸发 量 , 相 对 湿度 等 控制 水 分 条 件 的 因子 是 影响 
中 西部 地 区 NEP 的 主导 因子 ,特别 是 西部 地 区 降 
水 的 g 值 (gq=0.66) 最 大 。 而 对 于 东部 地 区 , 则 表现 出 
两 个 极为 明显 的 特征 ,一 方面 ,因子 数量 众多 且 整 
Aq 值 偏 小 ,经 筛选 后 的 驱动 因子 仍 有 11 个 ; 另 一 方 
面 ,因子 类 型 复杂 ,涉及 的 因子 包括 人 类 活动 .气象 
因素 ,以 及 地 形 因 素 。 这 表明 东部 地 区 NEP 与 其 驱 
动因 子 关系 十 分 复杂 ,难以 用 一 个 或 几 个 具有 代表 
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性 的 因子 说 明 NEP 的 驱动 力 ,但 综合 来 看 , 受 人 为 
干扰 的 影响 最 为 强烈 。 

2.4.2 交互 作用 探测 ”为 了 探究 不 同 驱动 因子 协同 
作用 下 对 NEP 分 布 的 影响 ,本 研究 利用 交互 作用 探 
测 器 分 析 不 同 驱动 因子 之 间 的 交互 作用 以 及 对 
NEP 空 间 分 异 的 影响 (图 6)。 结 果 表 明 ,(1) 任意 两 
个 驱动 因子 的 交互 作用 对 NEP 的 影响 均 大 于 单个 
子 的 独立 作用 ,表现 出 双 因 子 增强 或 非 线 性 增强 
效应 。 其 中 ,研究 区 内 ,温度 与 其 他 因子 的 交互 作 
用 全 部 表现 为 非 线 性 增强 (图 6d)。(2) 就 研究 区 而 
言 ,降水 同 其 他 各 驱动 因子 进行 交互 组 合 的 g 值 最 
大 ,表明 降水 是 影响 研究 区 NEP 的 主要 因素 。(3) 
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图 6 黄土 高 原 及 其 子 区 域 NEP 交互 作用 探测 结果 


Fig. 6 NEP interaction detection results of Loess Plateau and its sub-regions 
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中 西部 地 区 各 驱动 因子 的 交互 作用 g 值 普遍 较 高 ， 
以 降水 量 、 燕 发 量 等 代表 水 分 控制 条 件 的 因子 与 其 
他 因子 的 交互 作用 最 强 ; 而 东部 地 区 各 驱动 因子 的 
交互 作用 4 值 普遍 较 低 ,其 中 ,土地 利用 类 型 因子 与 


为 此 ,本 文通 过 与 其 他 学 者 的 NEP 佑 算 结果 进 
行 比 较 进 而 间接 评估 NEP 的 估算 精度 。 结 果 表 明 ， 
NEP 的 空间 分 布 及 数值 大 小 与 Chuai 等 , Zhang 
等 2 的 研究 结果 相似 ;NEP 趋 势 分 析 与 Ye 55^" BE 


a 


其 他 因子 的 交互 作用 最 强 。 
3 讨论 


3.1 研究 区 域 的 划分 

关于 黄土 高 原 整体 植被 恢复 已 开展 较 多 ,而 对 
于 高 原 内 部 的 差异 性 研究 相对 较 少 。 且 由 于 黄土 
高 原 自 然 地 理 环 境 的 复杂 性 ,学 者 们 往往 更 倾向 研 
究 地 理 区 划 引 起 的 差异 性 ,如 王 一 等 根据 《黄土 
高 原 地 区 综合 治理 规划 大 纲 (2010 一 2030 年 )》, 将 
黄土 高 原 按 地 貌 分 为 了 五 个 地 貌 区 域 ,进而 揭示 高 
原 内 部 植被 覆盖 度 的 差异 性 。 然 而 ,由 于 近 些 年 黄 
土 高 原 地 区 实施 了 一 系列 生态 恢复 措施 "i, 且 大 多 
以 行政 单元 为 主导 ,笔者 认为 侧重 于 人 类 活动 的 行 
政 区 划 间 的 差异 性 分 析 更 值得 关注 。 因 此 ,本 文 以 
大 的 省 级 行政 区 划 为 划分 单元 ,将 研究 区 按 经 度 方 
向 划分 为 西 . 中 、 东 三 个 子 区 域 ( 由 于 黄土 高 原 在 内 
蒙古 的 部 分 区 域 NPP 数 据 缺 失 ,本文 未 将 内 蒙古 纳 
和 信人 研究 范围 。 虽 未 完整 包涵 黄土 高 原 地 区 ,但 履 盖 
绝 大 部 分 地 区 ,并 突出 了 干 湿 变化 的 自然 差异 以 及 
由 行政 管理 导致 的 人 为 差异 ) ,在 一 定 程 度 上 既 能 
反映 出 政策 在 不 同行 政 区 实施 的 差异 性 ,又 尽 可 能 
保持 不 同 区 域 气候 等 自然 条 件 的 相似 性 ,从 而 更 能 
相对 完整 地 体现 生态 治理 工程 的 成 效 。 
3.2 _ NEP 估算 结果 的 不 确定 性 

全 球 或 区 域 尺度 的 NEP 佑 算 结果 的 精度 评价 
一 直 是 NEP 遥 感 估 算 的 难题 和 最 受 争议 的 环节 。 
稼 用 的 NEP 估 算 模 型 精度 验证 方法 通常 有 两 种 ,一 
种 是 实测 法 ,即将 估算 值 与 实际 测量 值 进行 对 比 ， 
以 验证 其 精度 ; 另 一 种 是 相对 法 ,即将 估算 的 空间 
分 布 图 与 其 他 学 者 的 研究 结果 或 冰 感 数据 产品 进 
行 对 比 ,以 评估 其 效果 。 尽 管 野外 调查 获得 的 实测 
数据 比较 可 靠 ,但 很 难 在 整个 研究 区 内 进行 大 范围 
比较 均匀 地 实地 调查 取样 。 此 外 ,估算 的 NEP 空 间 
分 辩 率 为 1 km, 时 间 跨 度 为 21 a, 生 物 量 分 布 的 时 空 
异 质 性 较 大 。 因 此 ,如 果 简 单 地 利用 有 限 的 实地 调 
查 数据 来 推算 整个 区 域 的 生物 量 , 则 可 能 产生 较 大 


究 结 果 相 符 ; 基 于 NEP 的 碳 源 / 汇 特征 区 域 的 分 布 及 
转换 格局 与 Zhang 等 2 的 结果 高 度 一 致 。 因 此 ,可 
以 认为 本 文 NEP 估 算 结 果 具 备 较 高 的 可 信和 度 。 
3.3 NEP 的 时 空 变化 

在 研究 时 段 内 ,黄土 高 原 NEP 呈 上 升 趋势 , 空 
间 上 ,研究 区 生态 环境 自 西北 向 东南 逐步 改善 ,这 
Ej Sun 等 ”的 人 研究 结果 相 一 至。 从 研究 区 内 部 看 ， 
中 部 地 区 是 区 内 NEP 变化 速率 最 大 的 区 域 ,其 NEP 
上 升 最 明显 的 区 域 位 于 关中 平原 和 陕 北 高 原 之 间 ， 
生态 治理 之 前 ,该 区 域 耕地 类 型 以 坡 耕 地 为 主 ， 
是 我 国 水 土 流 失 最 严重 的 地 区 之 一 所 。 由 于 退耕 
还 林 ( 草 ) 工 程 的 实施 ,该 区 域 碳 储量 和 生境 质量 逐 
步 提升 淖 。 而 榆林 地 区 北部 虽然 2020 年 NEP 相 较 
于 2000 年 已 经 有 所 增加 ,但 是 仍 为 碳 源 ,原因 可 能 
是 地 处 陕 北 高 原 ,靠近 毛乌素 沙 地 ,气候 干旱 ,土壤 
条 件 差 ,不 利于 植被 恢复 。 东 部 地 区 NEP 上 升 速率 
也 较 大 ,主要 分 布 在 水 土 流 失 严 重 的 山地 丘陵 区 ， 
且 该 区 域 气候 相对 湿润 ,利于 植被 恢复 和 生态 建 
设 。 而 西部 地 区 NEP 上 升 速率 相对 较 小 ,主要 分 布 
在 甘肃 和 陕西 的 交界 处 ,这 是 由 于 该 区 域内 大 多 处 
于 陇 中 高 原 , 且 降水 较 少 ,生态 治理 难度 大 ,而 靠近 
陕西 的 区 域 气候 条 件 相 对 于 整个 西部 地 区 而 言 较 
好 ,地 形 以 丘陵 为 主 呈 ,因此 ,NEP 上 升 较为 明显 。 
综 上 所 述 , 近 20 a, 通 过 坡 耕 地 退耕 还 林 还 草 等 人 为 
措施 ,黄土 高 原 植被 覆盖 明显 增加 呈 ,生态 碳 汇 功 
能 显著 增强 。 
3.4 NEP 的 驱动 因子 分 析 

研究 发 现 ,区 内 NEP 的 空间 分 异 特征 是 自然 环 
境 因 素 和 人 类 活动 因素 共同 作用 的 结果 。 已 有 研 
究 分 别 在 不 同时 空 尺 度 上 验证 了 气温 和 降水 皆 对 
NEP 变化 存在 影响 下 。 这 可 以 归 因 于 以 下 两 点 : 
(1) 气温 和 降水 是 估算 Rh 的 关键 因素 ,通过 影响 及 
间接 影响 NEP 的 大 小 ,(2) 气 温和 降水 代表 了 水 热 
条 件 , 通 过 影响 植物 生长 进而 影响 NEP 的 变化 。 通 
过 对 研究 区 2000 一 2020 年 气温 和 降水 年 际 变化 的 
研究 ,结果 表明 黄土 高 原 年 平均 气温 和 年 平均 降水 
量 分 别 增加 了 0.42 "和 91 mm( 图 7), 植 被 生长 条 件 
更 加 适宜 。 先 前 有 研究 表明 ,在 干旱 半 十 旱地 区 ， 
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(b) 降水 量 
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图 7 2000 一 2020 年 黄土 高 原 气 温和 降水 量 趋势 
Fig.7 Temperature and precipitation trends in the Loess Plateau from 2000 to 2020 


水 分 条 件 往往 是 影响 植被 恢复 的 重要 因素 %。 这 
一 发 现 与 本 文 结 果 一 致 ,文中 因子 探测 结果 也 显 
示 ,降水 .蒸发 量 以 及 相对 湿度 等 水 分 控制 因子 是 
区 内 NEP 分 布 的 主要 自然 驱动 因子 ,进一步 证 实 研 
究 区 NEP 分 布 主要 受 水 分 控制 。 

关于 人 类 活动 ,土地 利用 类 型 是 影响 黄土 高 原 
NEP 分 布 的 重要 驱动 因素 。 以 往 的 研究 表明 ,人 类 
活动 驱动 的 土地 利用 和 覆盖 变化 已 经 深刻 改变 了 
自然 生态 系统 的 结构 和 功能 ,不 仅 涉及 景观 ,还 改 
变 了 地 表 生 态 系 统 的 碳 通 量 ”。 土 地 利用 类 型 对 
NEP 的 影响 可 以 从 两 个 方面 来 解释 : (1) 不 同 土地 
利用 类 型 的 自然 属性 和 生态 结构 存在 差异 ,林地 在 
不 同 土地 利用 类 型 中 具有 最 强 的 碳 汇 能 力 ”, 这 与 
本 文中 植被 覆盖 度 ( 图 8) 以 及 NEP 变化 相似 ,(2) 
不 同 土地 利用 类 型 之 间 的 转移 必然 会 引起 NEP 的 
变化 ,因为 ,不 同 土地 利用 类 型 之 间 的 NEP 存 在 显 
著 的 差别 。 


(a) 2000 年 


人 


es 西部 地 区 
$ 中 部 地 区 
f ”东部 地 区 


交互 作用 探测 结果 表明 研究 区 NEP 的 空间 分 
布 受 降水 量 等 水 分 控制 因子 影响 明显 ,各 驱动 因子 
交互 组 合 的 解释 力 显著 高 于 单个 驱动 因子 。 同 时 ， 
本 人 研究 还 发 现 各 驱动 因子 与 温度 的 交互 作用 普遍 
呈现 非 线 性 增强 的 特征 (图 6d) ,对 此 ,笔者 目前 尚 
无 法 合理 解释 ,但 希望 这 一 现象 能 引起 相关 领域 学 
者 的 关注 。 此 外 ,在 研究 范围 内 ,三 个 子 区 域 的 
NEP 驱动 因子 存在 显著 的 空间 异 质 性 。 在 中 西部 
地 区 ,由 于 人 为 活动 干扰 相对 较 少 ,降水 量 、 共 发 量 
等 水 分 控制 因素 与 其 他 因素 的 交互 作用 决定 了 
NEP 的 空间 格局 。 而 在 东部 地 区 ,人 类 活动 对 区 内 
NEP 的 影响 程度 较 高 ,土地 利用 类 型 等 人 为 活动 因 
素 与 其 他 因素 的 交互 组 合 能 够 更 好 地 解释 NEP 的 
分 布 情况 。 

3.5 局 限 性 和 前 景 

植被 碳 固定 .土壤 呼吸 等 是 极其 复杂 的 生理 E 

态 过 程 , 通 过 模型 估算 均 存在 一 定 的 不 确定 性 ? , 仅 
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图 8 黄土 高 原 2000 年 和 2020 年 FVC 的 空间 格局 


Fig. 8 Spatial pattern of FVC in the Loess Plateau in 2000 and 2020 
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仅 基 于 NPP 的 碳 汇 能 力 评估 是 有 限 的 。 此 外 , 固 碳 
能 力 的 强 弱 是 一 个 长 期 动态 的 过 程 ,一 般 考 察 的 时 
间 序 列 越 长 ,评估 结果 越 可 靠 。 受 于 数据 来 源 的 限 
Hil ,本 研究 仅 选 取 2000—2020 4 NEP 作为 研究 对 
象 ,这 一 定 程 度 上 限制 了 NEP 及 其 驱动 因子 之 间 关 
系 的 深入 分 析 。 因 此 ,作者 下 一 步 计划 对 固 碳 能 
的 评估 基于 更 长 时 间 尺 度 .更 细 的 时 间 分 辩 率 (如 
月 和 季度 ) 以 及 更 高 精度 的 NEP 佑 算 模 型 ,从 而 考 
察 不 同 驱动 因子 下 回 碳 能 力 的 变化 ,为 区 域 生态 环 
境 人 恢复 提供 更 精确 、 细 致 的 科学 文 撑 。 
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4 结论 


-H 

本 研究 基于 NEP 佑 算 模 型 和 遥感 数据 ,进行 了 
2000 一 2020 年 黄土 高 原 NEP 的 时 空格 局 与 驱动 力 
分 析 , 主 要 结论 如 下 : 

(1) 2000 一 2020 年 ,人 研究 区 NEP 呈 东南 高 西北 
低 , 由 西北 向 东南 逐步 增加 趋势 ;不 同 区 域 之 间 差 
异 显著 ,NEP 上升 最 明显 的 区 域 多 分 布 在 山地 丘陵 
区 ,以 陕 . 甘 交界 处 最 为 突出 。 区 内 NEP 碳 汇 功能 
逐步 增强 , 近 20a 有 近 一 半 区 域 从 碳 源 特征 转换 为 
碳 汇 特征 区 域 。 

(2) 降水 是 影响 研究 区 NEP 空间 分 布 格局 的 主 
导 自 然 因素 ,改善 黄土 高 原水 分 条 件 是 进一步 改善 
黄土 高 原 NEP 的 关键 。 土 地 利用 类 型 是 影响 区 内 
NEP 空间 分 异 格 局 的 主导 人 为 因子 ,与 单个 因子 相 
比 ,自然 因子 和 人 为 因子 的 交互 作用 普遍 增强 了 对 
NEP 空 间 分 布 格局 的 影响 强度 。 

(3) 在 研究 区 内 ,三 个 子 区 域 NEP 的 驱动 因子 
具有 和 较 强 的 空间 分 异 特征 ,在 中 西部 地 区 ,NEP 受 
自然 环境 因素 影响 较 大 ,以 气候 因素 为 主 。 而 在 东 
部 地 区 ,NEP 受 气候 因素 ,地 形 地 貌 以 及 人 类 活动 
干扰 的 综合 影响 ,以 人 类 活动 干扰 为 主 。 

总 之 ,黄土 高 原 随 着 植被 的 恢复 ,具备 碳 汇 功 
能 的 区 域 面 积 显 著 增加 ,但 同时 也 存在 明显 的 空间 
STE. NEP 的 空间 分 布 格局 受众 多 因素 影响 ,其 
中 ,降水 和 土地 利用 类 型 的 影响 尤为 重要 。 
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Analysis of spatial and temporal patterns and drivers of local regional NEP in 
the Loess Plateau from 2000 to 2020 


PEI Hongze, ZHAO Yachao, ZHANG Tinglong 


(College of Resources and Environment, Northwest A & F University, Yangling 712100, Shaanxi, China) 


Abstract: The Loess Plateau is an important ecological barrier in China. In the context of global change, the 
carbon balance of the region is of great concern. In this study, based on MOD17A3HGF data, we analyzed the 
carbon source/sink characteristics of the local area of the Loess Plateau from 2000 to 2020 using the geostatistical 
model of soil respiration (GSMSR) model, trend analysis, difference analysis, and Geodetector, to reveal the 
spatial and temporal patterns of the region' s net ecosystem productivity (NEP) and its driving factors from 2000 
to 2020. Meanwhile, the study area was divided longitudinally into three subregions, west, center, and east, to 
compare the variability of drivers within different regions. The results showed that: (1) 49.69% of the localized 
area of the Loess Plateau shifted from a carbon source to a carbon sink within 21 years. NEP fluctuated upward 
with time, and was higher in the southeast than in the northwest, with a multi-year average of 12 g C-m^-a''* (2) 
Moisture condition was the main natural factor affecting the spatial distribution of NEP, while land use type was 
the main anthropogenic factor influencing the spatial distribution of NEP. In addition, the influence of the 
interaction between different factors on NEP was generally greater than that of single factors. The effects of 
interactions between different factors on NEP were generally greater than those of individual factors. (3) The 
factors driving NEP in the three subregions of west, central, and east have clear spatial differentiation 
characteristics, the central and western regions are more affected by climate, with precipitation, humidity, and 
other moisture conditions as the main ones. The eastern region is affected by a combination of factors such as 
topography, climate, and human activities, and anthropogenic interference, represented by land use type, is the 
strongest among these factors. This knowledge can serve as a basis for effective carbon management strategies 
and ecosystem conservation efforts in the region. 
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